
FAST 3-D SURFACE MULTIPLE PREDICTION 
CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims benefit of United States provisional patent application 
serial number 60/560,223, filed April 7, 2004, which is herein incorporated by reference. 

5 BACKGROUND OF THE INVENTION 
Field of the Invention 

[0002] Embodiments of the present invention generally relate to marine seismic 
surveying and, more particularly, to a method for attenuating the effect of surface 
multiples in seismic signals. 

1 0 Description of Related Art 

[0003] Seismic surveying is a method for determining the structure of subterranean 
formations in the earth. Seismic surveying typically utilizes seismic energy sources 
which generate seismic waves and seismic receivers which detect seismic waves. The 
seismic waves propagate into the formations in the earth, where a portion of the waves 
reflects from interfaces between subterranean formations. The amplitude and polarity 
of the reflected waves are determined by the differences in acoustic impedance 
between the rock layers comprising the subterranean formations. The acoustic 
impedance of a rock layer is the product of the acoustic propagation velocity within the 
layer and the density of the layer. The seismic receivers detect the reflected seismic 
waves and convert the reflected waves into representative electrical signals. The 
signals are typically transmitted by electrical, optical, radio or other means to devices 
which record the signals. Through analysis of the recorded signals (or traces), the 
shape, position and composition of the subterranean formations can be determined. 

[0004] Marine seismic surveying is a method for determining the structure of 
25 subterranean formations underlying bodies of water. Marine seismic surveying typically 
utilizes seismic energy sources and seismic receivers located in the water which are 
either towed behind a vessel or positioned on the water bottom from a vessel. The 
energy source is typically an explosive device or compressed air system which 
generates seismic energy, which then propagates as seismic waves through the body 
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of water and into the earth formations below the bottom of the water. As the seismic 
waves strike interfaces between subterranean formations, a portion of the seismic 
waves reflects back through the earth and water to the seismic receivers, to be 
detected, transmitted, and recorded. The seismic receivers typically used in marine 
5 seismic surveying are pressure sensors, such as hydrophones. Additionally, though, 
motion sensors, such as accelerometers may be used. Both the sources and receivers 
may be strategically repositioned to cover the survey area. 

[0005] Seismic waves, however, reflect from interfaces other than just those 
between subterranean formations, as would be desired. Seismic waves also reflect 

10 from the water bottom and the water surface, and the resulting reflected waves 
themselves continue to reflect. Waves which reflect multiple times are called 
"multiples". Waves which reflect multiple times in the water layer between the water 
surface above and the water bottom below are called "water-bottom multiples". Water- 
bottom multiples have long been recognized as a problem in marine seismic processing 

15 and interpretation, so multiple attenuation methods based on the wave equation have 
been developed to handle water-bottom multiples. However, a larger set of multiples 
containing water-bottom multiples as a subset can be defined. The larger set includes 
multiples with upward reflections from interfaces between subterranean formations in 
addition to upward reflections from the water bottom. The multiples in the larger set 

20 have in common their downward reflections at the water surface and thus are called 
"surface multiples". Figure 1 , discussed below, provides examples of different types of 
reflections. 

[0006] Figure 1 shows a diagrammatic view of marine seismic surveying. The 
procedure is designated generally as 100. Subterranean formations to be explored, 

25 such as 102 and 104, tie below a body of water 106. Seismic energy sources 108 and 
seismic receivers 1 10 are positioned in the body of water 106, typically by one or more 
seismic vessels (not shown). A seismic source 108, such as an air gun, creates 
seismic waves in the body of water 106 and a portion of the seismic waves travels 
downward through the water toward the subterranean formations 102 and 104 beneath 

30 the body of water 1 06. When the seismic waves reach a seismic reflector, a portion of 
the seismic waves reflects upward and a portion of the seismic waves continues 
downward. The seismic reflector can be the water bottom 1 12 or one of the interfaces 
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between two subterranean formations, such as interface 114 between fornriations 102 
and 104. When the reflected waves traveling upward reach the water/air interface at 
the water surface 116, a majority portion of the waves reflects downward again. 
Continuing in this fashion, seismic waves can reflect multiple times between upward 
5 reflectors, such as the water bottom 112 or formation interfaces below, and the 
downward reflector at the water surface 116 above, as described more fully below. 
Each time the reflected waves propagate past the position of a seismic receiver 110, 
the receiver 110 senses the reflected waves and generates representative signals. 

[0007] Primary reflections are those seismic waves which have reflected only once, 
10 from the water bottom 112 or an interface between subterranean formations, before 
being detected by a seismic receiver 110. An example of a primary reflection is shown 
in FIG. 1 by raypaths 120 and 122. Primary reflections contain the desired information 
about the subterranean formations which is the goal of marine seismic surveying. 
Surface multiples are those waves which have reflected multiple times between the 
15 water surface 116 and any upward reflectors, such as the water bottom 112 or 
formation interfaces, before being sensed by a receiver 110. An example of a surface 
multiple which is specifically a water bottom multiple is shown by raypaths 130, 132, 
134 and 136. The point on the water surface 116 at which the wave is reflected 
downward is generally referred to as the downward reflection point 133. The surface 
20 multiple starting at raypath 130 is a multiple of order one, since the multiple contains 
one reflection from the water surface 116. Two examples of general surface multiples 
with upward reflections from both the water bottom 112 and formation interfaces are 
shown by raypaths 140, 142, 144, 146, 148 and 150 and by raypaths 160, 162, 164, 
166, 168 and 170, Both of these latter two examples of surface multiples are multiples 
25 of order two, since the multiples contain two reflections from the water surface 116. In 
general, a surface multiple is of order i if the multiple contains i reflections from the 
water surface 116. Surface multiples are extraneous noise which obscures the desired 
primary reflection signal. 

[0008] Surface multiple attenuation is a prestack inversion of a recorded wavefield 
30 which removes all orders of all surface multiples present within the marine seismic 
signal. Unlike some wave-equation-based multiple-attenuation algorithms, surface 
multiple attenuation does not require any modeling of or assumptions regarding the 
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positions, shapes and reflection coefficients of the multiple-causing reflectors. Instead, 
surface multiple attenuation relies on the internal physical consistency between primary 
and multiple events that must exist in any properly recorded marine data set. The 
information needed for the surface multiple attenuation process is already contained 
5 within the seismic data. 

[0009] Various prior art methods have been tried for removal of surface multiples 
from recorded traces. It has been noted, for example, that the travel time for a water 
bottom multiple is a function of the "offset", the distance between the source and 
receiver, and the number of times the multiple reflects from the surface. For example, if 
10 the multiple reflects from the surface once before being received by the microphone 
and the offset is zero, the multiple's travel time is exactly twice that of the principal 
waves. This fact has been used in various schemes to remove multiples. 

[0010] other methods involve complex ray tracing schemes which generate a 
synthetic multiple wave and subtract it from the actual wave to obtain a supposedly 
15 multiple free record. However, these methods are very awkward in that they require 
significant knowledge of the subsea structure as well as the ocean bottom configuration 
before the synthetic wave can be generated. Similar synthetic multiples can be 
generated using more accurate methods not directly involving ray tracing, e.g., field 
propagation techniques, but again these require detailed knowledge of at least the 
20 ocean bottom, as well as the shape of the subsea interfaces, and so are not as 
practical as would be desired. 

[0011] Therefore, a need exists in the art for an improved method for removing the 
record of multiple surface reflection events from seismic records for seismic data 
processing purposes. 

25 SUMMARY OF THE INVENTION 

[0012] Embodiments of the present invention are generally directed to a method for 
predicting a plurality of surface multiples for a plurality of traces in a record of seismic 
data. In one embodiment, the method includes providing a plurality of target traces at a 
nominal offset and a nominal azimuth; selecting a plurality of pairs of input traces, 
30 wherein the midpoints of the input traces in each pair are separated by half the nominal 
offset and the azimuth of a line connecting the midpoints of the input traces in each pair 

4 



Attorney Docket No.: WGEOBRsPCT 



is equal to the nominal azimuth; convolving the selected pairs of input traces to 
generate a plurality of convolutions; and applying a three dimensional operator to the 
convolutions. 

[0013] In another embodiment, the method includes dividing a plurality of target 
5 traces into one or more groups according to offsets; dividing each group into one or 
more subgroups according to azimuths; selecting a first subgroup having a first nominal 
offset and a first nominal azimuth; selecting a plurality of pairs of input traces, wherein 
the midpoints of the input traces in each pair are separated by half the first nominal 
offset and the azimuth of a line connecting the midpoints of the input traces in each pair 
10 is equal to the first nominal azimuth; convolving the selected pairs of input traces to 
generate a plurality of convolutions; and applying a three dimensional operator to the 
convolutions. 

[0014] In yet another embodiment, the method includes providing a plurality of target 
traces at a nominal offset; predicting a plurality of two dimensional surface multiples for 
15 a plurality of input subsurface lines; and applying a two dimensional operator to the 
predicted two dimensional surface multiples. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] So that the manner in which the above recited features of the present 
invention can be understood in detail, a more particular description of the invention, 
20 briefly summarized above, may be had by reference to embodiments, some of which 
are illustrated in the appended drawings. It is to be noted, however, that the appended 
drawings illustrate only typical embodiments of this invention and are therefore not to 
be considered limiting of its scope, for the invention may admit to other equally effective 
embodiments. 

25 [0016] Figure 1 illustrates a diagrammatic view of marine seismic surveying. 

[0017] Figure 2 a flow diagram of a method for performing a three dimensional 
surface multiple prediction in accordance with one or more embodiments of the 
invention. 
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[0018] Figure 3 illustrates a flow diagram of a method for performing a three 
dimensional surface multiple prediction in accordance with one or more embodiments 
of the invention, 

[0019] Figure 4 illustrates a plan view of an acquisition geometry in accordance with 
5 one or more embodiments of the invention. 

[0020] Figure 5 illustrates a computer network into which various embodiments of 
the invention may be implemented. 

DETAILED DESCRIPTION 

[0021] Figure 2 illustrates a flow diagram of a method 200 for performing a three 
10 dimensional surface multiple prediction in accordance with one or more embodiments 
of the invention. At step 210, a set of target traces is divided into one or more groups 
according to offsets. Each group contains target traces with offsets within a specified 
range, where an offset is defined as the horizontal distance between a source and a 
receiver. The set of target traces defines the locations at which the multiples are to be 
15 predicted. The offset ranges span from the shortest offset in the set of target traces, 
which typically corresponds to the receiver located closest to the back of the vessel, 
e.g., about 100 meters, to the longest offset in the set of target traces, which typically 
corresponds to the receiver furthest from the back of the vessel, e.g., about 6000 
meters. Each offset range has a nominal offset value, which is typically the central 
20 offset value in that range. For example, target traces with offsets between 100 meters 
and 150 meters may be organized into one range with a nominal offset value of 125 
meters, and target traces with offsets between 150 meters and 200 meters may be 
organized into another range with a nominal offset value of 1 75 meters. 

[0022] At step 220, a group, such as a first group, of the target traces with offsets 
25 within a first range of offsets is selected. At step 230, the selected group of the target 
traces is divided into one or more subgroups of target traces according to azimuths. An 
azimuth is defined as the angle between the line that connects the source and the 
receiver and some fixed direction, which is typically the in-line direction. Each azimuth 
range has a nominal value, which is typically the central value of that range. In this 
30 manner, all the target traces within a subgroup are at the same or similar azimuths (or 



6 



Attorney Docket No.: WGEflBPsPCT 



have the same nominal azimuth) and the same or similar offsets (or have the same 
nominal offset). At step 240, a subgroup, such as a first subgroup, of target traces is 
selected. The first subgroup of target traces is at the first nominal offset and the first 
nominal azimuth. 

5 [0023] At step 250, a set of input traces is selected. The input traces are the traces 
from which the predicted multiples are computed. In one embodiment, a set of input 
traces at half of the nominal offset for the selected subgroup of target traces is 
selected. However, a set of input traces at any offset range may be selected. In 
another embodiment, the input traces are interpolated and regularized. Alternatively, 
10 several offsets may be merged and/or sub-stacked to improve sampling and signal to 
noise ratio. 

[0024] At Step 260, the selected input traces are preconditioned to simulate zero 
offset traces. The selected input traces may be preconditioned using many techniques 
well known in the art, such as a moveout correction, a full migration/demigration and 
15 the like. The input traces may be located by their midpoints. In this manner, each 
preconditioned input trace simulates a trace that would have been recorded with a 
source and a receiver at the midpoint location. 

[0025] At step 270, each pair of selected preconditioned input traces is convolved to 
generate a set of convolutions. The pairs are selected such that the separation 
between the midpoints of the traces in each pair is equal to half the nominal offset for 
the selected subgroup of target traces, and the azimuth of the line connecting the 
midpoints of the traces in each pair is equal to the nominal azimuth of the selected 
subgroup of target traces. Each convolution may be located at the midpoint of the line 
connecting the midpoints of the two traces in each pair. The input traces may also be 
interpolated to form the convolutions. Actual or possible structural dips may also be 
taken into account in forming the convolutions. 

[0026] At Step 280, a three dimensional operator is applied to the set of convolutions 
and the result is located at the midpoints of the selected subgroup of target traces. In 
one embodiment, the three dimensional operator is a three dimensional demigration 
30 operator with an offset equal to half of the nominal offset of the selected subgroup of 
target traces and a velocity equal to half the water velocity. In such a case, the result 
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may be referred to as demigrated convolutions. Another example of the three 
dimensional operator is a poststack (zero-offset) demigration followed by an inverse 
Dip Moveout (DMO) and inverse moveout correction. In this manner, the surface 
multiples for the selected subgroup of target traces may be predicted. In one 
embodiment, the result may be corrected from the nominal offset and azimuth to the 
actual offset and azimuth. For instance, the offset correction may be performed using a 
differential moveout correction, and the azimuth correction may be performed via 
interpolation from adjacent azimuths. 

[0027] Alternatively, the Input traces may be convolved at their original offset and 
the convolutions may be corrected prior to applying the three dimensional operator. The 
correction may be performed using a moveout correction with half the velocity (or twice 
the offset). 

[0028] At step 285, a determination is made as to whether another subgroup, such 
as a second subgroup, of target traces with a second nominal azimuth exists. If the 
answer is in the affirmative, then processing returns to step 240, at which the second 
subgroup of target traces is selected. Processing then continues through steps 250- 
280 at which a set of input traces (which may or may not be the same as the first set) is 
processed using the nominal azimuth of the second subgroup of target traces in order 
to predict multiples for those traces. Steps 240-285 continue until all the subgroups 
within the selected group have been processed. 

[0029] If the answer is in the negative, then processing continues to step 290, at 
which a determination is made as to whether another group, such as a second group, 
of target traces with a second nominal offset exists. If the answer is in the affirmative, 
then processing returns to step 220, at which the second group of target traces with the 
second nominal offset is selected. Processing then continues through steps 230-285 at 
which one or more sets of input traces are processed using the nominal offset for the 
second group of target traces in order to predict multiples for those traces. Steps 220- 
290 continue until all the groups within the set of target traces have been processed. 

[0030] In one embodiment, method 200 may be used to predict the three 
dimensional surface multiples when substantial streamer (cable) feathering exists. In 
another embodiment, method 200 may be used to predict multiples at a specified 
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azimuth and may therefore be made to account for azimuth variations in the seismic 
data. These azimuth variations can be a significant source of error in predicted 
multiples. As such, in accordance with another embodiment of the invention, method 
200 may be used to predict errors in a manner similar to two dimensional and three 
5 dimensional surface multiple prediction methods described in US Patent Application 
Serial No. 1 0/668,927, which is incorporated herein by reference. 

[0031] Figure 3 illustrates a flow diagram of a method 300 for performing a three 
dimensional surface multiple prediction in accordance with one or more embodiments 
of the invention. At step 310, a set of target traces is divided into one or more groups 
10 of traces according to offsets. The set of target traces defines the locations at which 
the multiples are to be predicted. Each group contains target traces with offsets within 
a specified range. As with the target traces described with reference to step 210, each 
offset range may have a nominal offset value, which is typically the central offset value 
in that range. 

[0032] At step 320, a record of input seismic traces is preconditioned. That is, the 
record of input seismic traces is separated into input subsurface lines (SSLs) and each 
input SSL is regularized according to conventional regularization methods known by 
persons of ordinary skill in the art. Once regularized, the crossline offset between each 
source and receiver is zero and the inline offset between each source and receiver is 
regular. Once the record of input seismic traces has been regularized to input SSL's, 
the regularized traces are extrapolated according to conventional extrapolation 
methods known by persons of ordinary skill in the art. Once extrapolated, the gap 
between each source and the receiver nearest to the source on each input SSL is filled 
with extrapolated receivers. As a result, each input SSL has traces with zero crossline 
offset and regularly increasing inline offsets starting from zero. 

[0033] At step 330, the two dimensional surface multiples for each input SSL are 
predicted using a conventional two dimensional prediction algorithm known by persons 
of ordinary skill in the art. At step 340, the predicted two dimensional surface multiples 
are sorted into one or more planes (or groups) according to offsets. Each plane 
30 contains multiples within a specified range of offsets. Each offset range may have a 
nominal offset value, which is typically the central offset value in that range. At step 
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345, a plane of the predicted two dimensional surface multiples with a nominal offset, 
such as a first nominal offset, is selected. 

[0034] At step 350, a two dimensional operator is applied to the plane of predicted 
two dimensional surface multiples and the result is located at the midpoints of a group 
5 of target traces with the first nominal offset. The two dimensional operator may be a 
two dimensional demigration operator with a velocity equal to half the water velocity or 
half the primary velocity. In this manner, the two dimensional multiple prediction takes 
into account variations in the input SSL's in the in line direction and the two dimensional 
operator takes into account variations in the input SSL's in the cross line direction. 

10 [0035] At Step 360, a determination is made as to whether a plane of predicted two 
dimensional surface multiples at a second nominal offset exists. If the answer is in the 
affirmative, then processing returns to step 345, at which a plane of predicted two 
dimensional surface multiples at a second nominal offset is selected. Processing then 
continues to step 350, at which a two dimensional operator is applied to the plane of 

15 predicted two dimensional surface multiples at the second nominal offset and the result 
is located at the midpoints of a group of target traces at the second nominal offset. 
Steps 345-360 continue until all of the predicted two dimensional surface multiples at 
various nominal offsets have been processed. In one embodiment, method 300 may 
be used to predict the three dimensional surface multiples when the streamers (or 
20 cables) are all substantially parallel with each other with negligible feathering. 

[0036] The following paragraphs provide mathematical derivations for methods 200 
and 300 in accordance with one or more embodiments of the invention. Figure 4 
illustrates a plan view of an acquisition geometry in accordance with one or more 
embodiments of the invention. Surface multiples for a trace (S, R), with source at S 
25 and receiver at R, are to be predicted. M and h are defined as the midpoint and offset 
of (S, R) respectively. X is defined as a potential downward reflection point (DRP) for 
the surface multiples. 

[0037] A three dimensional surface multiple prediction may be performed by 
convolving trace (S, X) with trace (X, R), and summing these convolutions over all 
30 possible X. In order to do this, traces (S, X) and (X, R) generally need to be estimated 
from recorded traces. One way to estimate trace (S, X) is to apply a differential 
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moveout correction to a trace with a similar midpoint, Ms, from the recorded dataset 
and with a similar offset. Trace (X, R) may be estimated in like manner. The 
differential moveout applied to the trace used to estimate (S, X) depends not only on its 
offset, but also on the offset from S to X, denoted xs in Figure 4, Typically, a given 
5 recorded trace will only be used to estimate one trace for each subsurface line (SSL) in 
the prediction process, but it will be used on every SSL within a given aperture of the 
recorded trace. As such, the differential moveout correction and convolution may be 
repeated for every SSL for which the prediction is required. 

[0038] The order of the moveout correction and the convolution may be swapped, 
10 i.e. traces with midpoints at Ms and Mr are convolved before the moveout correction is 
applied, and the resultant error may be corrected after the convolution. The recorded 
traces may be processed to approximate (kinematically) zero-offset traces at the 
midpoint locations. For a simple, first-order multiple from a horizontal reflector at time 
in a constant velocity medium with velocity the primary reflections occur at time to on 
15 the zero-offset traces at Ms and Mr, and at time 2<b on their convolution. The travel 
time for the primary reflection on trace (S, X) is given by: 

t{S,X) = [to^^{xs/vff (1) 

and similarly for the trace (X, R), such that the travel time on the convolution of these 
traces is given by 

20 m = -f (xs/v)T + Ih^ + (xB/\dT (2) 

[0039] Accordingly, a post-convolution correction is performed to map an event at 
time 2to to time f(X). At this point, a constant velocity demigration operator with velocity 
l/maps a migrated time fm to a demigrated time id given by: 

fd = [(U2f + {Xs/V)T + [{U2f + {Xb/V)T (3) 

25 where Xs and Xr are the distances from the migrated location to the demigrated source 
and receiver locations, respectively. Equation (3) has a similar form to the post- 
convolution correction. Thus, if = 2fo and fd = ^(X), then equation (3) becomes 

t(X) = + {Xs/V)T + [to^ + (Xb/V)T (4) 
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[0040] Equation (4) would be identical to the post-convolution correction If the 
convolved trace is located at X, and Xs = xs, Xr = xr and V= v. However, the location 
X depends on S and R, and hence the demigration would be repeated whenever a 
change in (S, R) creates a change in X. The convolved trace may be located at the 
"migrated location" Mx, the midpoint of Ms and Mr, which does not depend on S or R, 
but only on Ms and Mr. If S' and R' are defined to be the midpoints of (S, M) and (M, 
R) respectively, then the distances from Mx to S' and R' are xs/2 and xr/2 respectively. 
Hence, demigration from Mx to (S',R') yields a demigrated traveltime given by equation 
(3) with Xs = xs/2 and Xr = xr/2, which provides the desired demigrated traveltime if 
velocity V= v/2, 

[0041] Notably, (Ms, Mr) is parallel to (S, R) and the offset from Ms to Mr is h/2. 
Accordingly, if the offset and azimuth of (S, R) are fixed, then the offset and azimuth of 
(Ms. Mr) are also known. If Mx is determined and a volume of data that has been 
kinematically mapped to zero offset is available, then Ms and Mr may be determined, 
and the corresponding traces may be convolved from this volume, thereby placing the 
convolution into a new volume at location Mx. Repeating this operation for all Mx, and 
then demigrating the volume at constant offset {h/2) and azimuth (as defined 
previously) with velocity v/2 yields the predicted multiples for the entire offset plane. 

[0042] One distinction between the embodiments of the invention and the prior art is 
that the traces are convolved at fixed (often zero) offset before the offset is corrected. 
The improved performance accomplished by various embodiments of the invention is 
attained by making some approximations in the derivation of the algorithm. Theory and 
tests have shown that, in many cases, the predicted multiples are still sufficiently 
accurate that they may be adaptively subtracted. When crossline dip effects are 
significant, multiples predicted using various embodiments of the invention will be much 
more accurate than those predicted by a two dimensional surface multiple algorithm. 

[0043] Figure 5 illustrates a computer network 500, into which embodiments of the 

invention may be implemented. The computer network 500 includes a system 

computer 530, which may be implemented as any conventional personal computer or 

workstation, such as a UNIX-based workstation. The system computer 530 is in 

communication with disk storage devices 529, 531, and 533, which may be external 

hard disk storage devices. It is contemplated that disk storage devices 529, 531 , and 
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533 are conventional hard disk drives, and as such, will be implemented by way of a 
local area network or by remote access. Of course, while disk storage devices 529, 
531, and 533 are illustrated as separate devices, a single disk storage device may be 
used to store any and all of the program instructions, measurement data, and results as 
desired. 

[0044] In one embodiment, seismic data from hydrophones are stored in disk 
storage device 531. The system computer 530 may retrieve the appropriate data from 
the disk storage device 531 to perform the 3-D surface multiple prediction according to 
program instructions that correspond to the methods described herein. The program 
instructions may be written in a computer programming language, such as C++, Java 
and the like. The program instructions may be stored in a computer-readable memory, 
such as program disk storage device 533. Of course, the memory medium storing the 
program instructions may be of any conventional type used for the storage of computer 
programs, including hard disk drives, floppy disks, CD-ROMs and other optical media, 
magnetic tape, and the like. 

[0045] According to the preferred embodiment of the invention, the system computer 

530 presents output primarily onto graphics display 527, or alternatively via printer 528. 
The system computer 530 may store the results of the methods described above on 
disk storage 529, for later use and further analysis. The keyboard 526 and the pointing 
device (e.g., a mouse, trackball, or the like) 525 may be provided with the system 
computer 530 to enable interactive operation. 

[0046] The system computer 530 may be located at a data center remote from the 
survey region. The system computer 530 is in communication with hydrophones (either 
directly or via a recording unit, not shown), to receive signals indicative of the reflected 
seismic energy. These signals, after conventional formatting and other initial 
processing, are stored by the system computer 530 as digital data in the disk storage 

531 for subsequent retrieval and processing in the manner described above. While 

Figure 5 illustrates the disk storage 531 as directly connected to the system computer 

530, it is also contemplated that the disk storage device 531 may be accessible through 

a local area network or by remote access. Furthermore, while disk storage devices 

529, 531 are illustrated as separate devices for storing input seismic data and analysis 

results, the disk storage devices 529, 531 may be implemented within a single disk 
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drive (either together with or separately from program disk storage device 533), or in 
any other conventional manner as will be fully understood by one of skill in the art 
having reference to this specification. 

[0047] While the foregoing is directed to embodiments of the present invention, 
5 other and further embodiments of the invention may be devised without departing from 
the basic scope thereof, and the scope thereof is determined by the claims that follow. 
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What Is Claimed Is: 

1. A method for predicting a plurality of surface multiples for a plurality of traces in 
a record of seismic data, comprising: 

5 providing a plurality of target traces at a nominal offset and a nominal azimuth; 

selecting a plurality of pairs of input traces, wherein the midpoints of the input 
traces in each pair are separated by half the nominal offset and the azimuth of a line 
connecting the midpoints of the input traces in each pair is equal to the nominal 
azimuth; 

10 convolving the selected pairs of input traces to generate a plurality of 

convolutions; and 

applying a three dimensional operator to the convolutions. 

2. The method of claim 1, further comprising preconditioning the input traces to 
simulate zero offset traces. 

15 3. The method of claim 1, further comprising preconditioning the input traces to 
simulate zero offset traces by a moveout correction. 

4. The method of claim 1 , wherein the target traces provide the locations at which 
the surface multiples are predicted. 

20 

5. The method of claim 1, wherein the nominal offset is a central offset value within 
a range of offsets. 

6. The method of claim 1 , wherein the nominal azimuth is a central azimuth value 
25 within a range of azimuths. 

7. The method of claim 1 , wherein selecting the pairs of input traces comprises 
selecting the input traces at half the nominal offset. 

30 8. The method of claim 1 , further comprising Interpolating and regularizing the input 
traces. 
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9. The method of claim 1 , further comprising locating each convolution at the 
midpoint of the line connecting the midpoints of the input traces in each pair. 

10. The method of claim 1 , wherein applying the three dimensional operator to the 
convolutions comprises locating the result of the three dimensional operator application 
at the midpoints of the target traces. 

1 1 . The method of claim 1 0, further comprising correcting the result to an actual 
offset and an actual azimuth. 

12. The method of claim 1 , wherein the three dimensional operator is a three 
dimensional demigration operator having an offset equal to half the nominal offset and 
a velocity equal to one of half the water velocity or half of a multiple velocity function. 

13. A method for predicting a plurality of surface multiples for a plurality of traces in 
a record of seismic data, comprising: 

dividing a plurality of target traces into one or more groups according to offsets; 
dividing each group into one or more subgroups according to azimuths; 
selecting a first subgroup having a first nominal offset and a first nominal 
azimuth; 

selecting a plurality of pairs of input traces, wherein the midpoints of the input 
traces in each pair are separated by half the first nominal offset and the azimuth of a 
line connecting the midpoints of the input traces in each pair is equal to the first nominal 
azimuth; 

convolving the selected pairs of input traces to generate a plurality of 
convolutions; and 

applying a three dimensional operator to the convolutions. 

14. The method of claim 13, wherein the target traces within each group are at a 
nominal offset. 

15. The method of claim 14, wherein the nominal offset is a central offset value 
within a range of offsets. 
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1 6. The method of claim 1 3, wherein the target traces within each subgroup are at a 
nominal offset and a nominal azimuth. 

17. The method of claim 15, wherein the nominal azimuth is a central azimuth value 
within a range of azimuths. 

18. The method of claim 13, further comprising preconditioning the input traces to 
simulate zero offset traces. 

1 9. The method of claim 1 3, further comprising locating each convolution at the 
midpoint of the line connecting the midpoints of the input traces in each pair. 

20. The method of claim 1 3, wherein applying the three dimensional operator to the 
convolutions comprises locating the result of the three dimensional operator application 
at the midpoints of the target traces. 

21 . The method of claim 13, wherein the three dimensional operator is a three 
dimensional demigration operator having an offset equal to half the first nominal offset 
and a velocity equal to one of half the water velocity or half of a multiple velocity 
function. 

22. A method for predicting a plurality of surface multiples for a plurality of traces in 
a record of seismic data, comprising: 

providing a plurality of target traces at a nominal offset; 

predicting a plurality of two dimensional surface multiples for a plurality of input 
subsurface lines; and 

applying a two dimensional operator to the predicted two dimensional surface 
multiples. 

23. The method of claim 22, further comprising, prior to predicting the two 
dimensional surface multiples, preconditioning the input traces into the input subsurface 
lines. 
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24. The method of claim 22, further comprising sorting the predicted two 
dimensional surface multiples into one or more planes according to offsets, wherein 
each plane has an associated nominal offset, 

25. The method of claim 24, further comprising selecting a plane of predicted two 
dimensional surface multiples having a first nominal offset. 

26. The method of claim 22, wherein applying the two dimensional operator 
comprises locating the result of the two dimensional operator application at the 
midpoints of the target traces. 

27. The method of claim 22, wherein the two dimensional operator is a two 
dimensional demigration operator with a velocity equal to one of half the water velocity 
or half of a multiple velocity function. 
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ABSTRACT OF THE DISCLOSURE 

A method and apparatus for predicting a plurality of surface multiples for a 
plurality of traces in a record of seismic data. In one embodiment, the method includes 
providing a plurality of target traces at a nominal offset and a nominal azimuth; 
selecting a plurality of pairs of input traces, wherein the midpoints of the input traces in 
each pair are separated by half the nominal offset and the azimuth of a line connecting 
the midpoints of the input traces in each pair is equal to the nominal azimuth; 
convolving the selected pairs of input traces to generate a plurality of convolutions; and 
applying a three dimensional operator to the convolutions. 
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